Resonant electromagnetic emission from intrinsic Josephson-junction stacks with 
laterally modulated Josephson critical current 



A. E. Koshelev 

Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439 

L. N. Bulaevskii 
Los Alamos National Laboratory, Los Alamos, New Mexico 87545 
(Dated: February 3, 2008) 

Intrinsic Josephson-junction stacks realized in mesas fabricated out of high-temperature super- 
conductors may be used as sources of coherent electromagnetic radiation in the terahertz range. 
The major challenge is to synchronize Josephson oscillations in all junctions in the stack to get 
significant radiation out of the crystal edge parallel to the c axis. We suggest a simple way to solve 
this problem via artificially prepared lateral modulation of the Josephson critical current identical in 
all junctions. In such a stack phase oscillations excite the in-phase Fiske mode when the Josephson 
frequency matches the Fiske-resonance frequency which is set by the stack lateral size. The powerful 
almost standing electromagnetic wave is excited inside the crystal in the resonance. This wave is 
homogeneous across the layers meaning that the oscillations are synchronized in all junctions in the 
stack. We evaluate behavior of the I-V characteristics and radiated power near the resonance for 
arbitrary modulation and find exact solutions for several special cases corresponding to symmetric 
and asymmetric modulations of the critical current. 



I. INTRODUCTION 

Josephson junctions are natural voltage-to-frequency 
converters, since a finite voltage drop across the junc- 
tion always leads to oscillating current with frequency 
proportional to the voltage (ac Josephson effect 1 ). This 
fundamental property suggested an attractive possibility 
to use the ac Josephson effect for developing of voltage- 
tunable generators of electromagnetic waves. Radiation 
from a Josephson junction directly into the waveguide in 
the microwave frequency range has been indeed detected 
a long time ago 2 ^. However, the typical detected radi- 
ated power ~1 pW occurred to be too small for practical 
applications. 

A natural route to enhance this power is to use large 
arrays of junctions. This possibility has been extensi vely 
explored by several experimental groups, see reviews^. 
When all junctions oscillate in phase, the total emitted 
power is expected to be proportional to the square of the 
total number of junctions in the array. Inevitable varia- 
tions of junction parameters, however, may cause varia- 
tions of the oscillating frequencies leading to dcsynchro- 
nization and dramatic drop in emission power. There- 
fore, the major challenge is to synchronize oscillations in 
all junctional One way to solve this problem is to cou- 
ple the junctions with a resonant cavity. The efficient 
synchronization in such systems has been demonstrated 
experimentally*^ and has been extensively studied in sev- 
eral simulation papers^'. 

Demonstration of intrinsic Josephson effect in high- 
temperature superconductors^ opened a completely new 
route to developing electromagnetic sources. Lay- 
ered high-temperature superconducting materials, such 
as Bi 2 Sr 2 CaCu208 (BSCCO), are composed of super- 
conducting Cu02 layers coupled by Josephson inter- 
action. Intrinsic Josephson effect has been exten- 



sively investigated in the past decade and most "clas- 
sical" ac and dc Josephson phenomena have been ob- 
served, including Fraunhofcr magnetic oscillations of crit- 
ical curr ent in small-size samples^, Josephson plasma 
resonanc e 1 11 ^ 12 !, Shapiro steps in current-voltage charac- 
teristics induce d by ext ernal microwave irradiatiod 13 ! 14 !, 
Fiske resonance d 15 ! 16 ! 17 -!, etc. Therefore a small-size mesa 
fabricated out of this material represents a natural one- 
dimensional array of Josephson junctions. A large value 
of superconducting gap, up to 60 meV, allows to bring 
the Josephson frequency into the practically important 
terahertz range. Due to atomic nature of the intrinsic 
junctions, it is much easier to prepare large arrays of vir- 
tually identical junctions than in the case of artificially 
fabricated arrays. Nevertheless, the same challenge to 
synchronize oscillations in all junctions also remains for 
this system. 

Electromagnetic waves propagate inside large-size lay- 
ered superconductor in the form of plasma modes. In zero 
magnetic field the minimum frequency of these waves cor- 
responding to homogeneous oscillations is given by the 
Josephson plasma frequency. The in-plane velocity of 
the mode strongly depends on the c-axis wave vector, q z . 
The highest velocity corresponds to the in-phase mode, 
q z = 0, and the lowest velocity corresponds to the an- 
tiphase mode, q z = n/s. 

A stack of the intrinsic Josephson junctions with lat- 
eral size smaller than the decay length of electromagnetic 
wave behaves as a cavity, ft is characterized by a spec- 
trum of Fiske resonant modes corresponding to excitation 
of almost standing electromagnetic waves-^. Frequencies 
of these modes strongly depend on the wave vector per- 
pendicular to the layer direction, with the maximum fre- 
quency corresponding to the homogeneous in-phase mode 
in all junctions and the minimum frequency correspond- 
ing to the antiphase mode. At the stack edge the electro- 
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magnetic waves excited inside the intrinsic junctions are 
converted into electromagnetic waves propagating into 
dielectric media outside the crystal. To use the stack as 
a coherent source of such radiation, it would be desir- 
able to excite the in-phase resonance mode. Then the 
radiation power is proportional to the square of the total 
number of junctions positioned at distances smaller than 
the wavelength of radiation. However, to synchronize 
the intrinsic junctions one needs to have strong enough 
coupling between them. In the simplest case of a homo- 
geneous stack at zero magnetic field, the Josephson oscil- 
lations typically interact very weakly inside the crystal. 

A very popular way to achieve coherent Josephson os- 
cillations in the whole stack is to apply magnetic field 
along the layers. The magnetic field promotes strong 
inductive interaction between the neighboring junctions. 
Large magnetic field generates a dense Josephson vortex 
(JV) lattice homogeneously filling all junctions. In the 
case of large-size system, interaction between the static 
JV arrays in neighboring junctions leads to formation of 
the triangular lattice corresponding to the 7r phase shift 
between the phases in the neighboring junctions. The 
Fiske resonances excited by t he mov ing JV lattice have 
been observed experimental!}^ 16 * 17 ! . Due to its trian- 
gular ground state, the JV lattice typically excites the 
antiphase modes. Only a very weak outside radiation at 
double Josephson frequency is expected in this case.^ 

To achieve a powerful emission it would be desirable 
to prepare aligned rectangular arrangement of JVs. Such 
configuration is expected at certain conditions in small- 
size mesas due to interaction with edges. At small lat- 
tice velocities transition from triangular to rectangular 
configuration with increasing magnetic field has been ob- 
served experimentally as a crossover from 3>o/2- to <&o- 
periodic magnetic oscillations of the flux-flow voltage^. 
The transitions between static configurations have been 
studied theoretically in Ref. 21. The possibility to have 
the aligned configuration at large velocities is an open 
issue. Recent large-scale numerical simulations 22 -! suggest 
that excitation of the resonance in-phase mode and in- 
teraction via the radiation electromagnetic field may pro- 
mote alignment of JVs. 

A possibility to use a mesa with small lateral size and a 
very large number of junctions (about 10 ) in zero mag- 
netic field as a source of terahertz radiation has been pro- 
posed in Ref. [521 In such a design oscillations in different 
junctions are synchronized by the external electromag- 
netic radiation field generated by the oscillations them- 
selves. Small lateral size increases the strength of interac- 
tion due to the radiation field and allows to avoid exces- 
sive heating. Such a source is frequency-tunable with the 
maximum power conversion efficiency about 30%. The 
obvious technological challenge of this design is a require- 
ment to fabricate a mesa with such a large number of 
almost identical junctions. 

In this paper we explore a different way to excite res- 
onance mode and synchronize oscillations in a junction 
stack. We propose to use a stack with strongly modulated 



Josephson critical current (JCC) with modulation identi- 
cal in all junctions. Such a modulation dramatically en- 
hances coupling between the Josephson oscillations and 
the in-phase resonance modes. For a single junction such 
a mechanism of excitation of the Fiske resonances has 
been considered in Ref. EH 

The frequencies of in-phase resonance modes are set by 
the lateral size of the mesa, L, 
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where e c is the c-axis dielectric constant. In particular, 
assuming e c = 12, the resonance at uj\/2tt = 1 THz for 
the fundamental mode, m = 1, takes place for a mesa 
with width 43 /im. In the resonance a powerful almost 
standing wave is excited by Josephson oscillations which 
synchronizes oscillations in the whole stack. Such syn- 
chronization function of the cavity mode does not exist 
in a single junction!^ In such a design the frequency 
tunability in a single device is sacrificed in favor of larger 
power and better efficiency. We consider several spe- 
cific cases of modulation, both symmetric and asymmet- 
ric, which correspond to excitation of the mode with the 
wavelength equal to L and 2L. For simplicity, we assume 
that the mesa size along the field direction, L y , is larger 
than the wavelength of outcoming electromagnetic wave, 
A w (0.3 mm for 1 THz). The calculation can be straight- 
forwardly generalized to the opposite case L y < A w . We 
calculate the I-V dependences, radiated electromagnetic 
power, and power conversion efficiency for such systems. 

Recently, the resonant electromagnetic emission from 
the mesas fabricated out of underdoped BSCCO crystals 
has been detected experimentally^. The resonance fre- 
quencies vary from 0.4 to 0.85 THz, and they increase 
roughly inversely proportional to the mesa widths. The 
origin of the observed resonances is most probably due to 
the mechanism described in this paper, even though no 
JCC modulation has been introduced deliberately Sup- 
pression of superconductivity near the edges during the 
fabrication process occurs to be sufficient for excitation 
of the resonances. One can expect that deliberately in- 
troduced JCC modulation will significantly enhance the 
amplitude of the resonance and radiation power. 

The paper is organized as follows. In section [H] we 
present the equation and boundary conditions for the 
oscillating phase when it is identical in all junctions. 
Derivation of the boundary condition for this case is sum- 
marized in Appendix [A] We also present the radiation 
power in terms of the oscillating phase and power con- 
version efficiency. In section |III| we derive the energy 
balance relations in the dynamic state. In section IV 



using these relations, we analyze the behavior near the 
resonances and derive approximate results for resonant 
enhancements of the current, radiated power, and power 
conversion efficiency. Appendix [5] presents derivation of 
the radiation losses for the resonance mode in a thin rect- 
angular mesa. In section [V] we consider several special 
cases of modulation for which the problem allows for ex- 
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act analytical solutions, see Fig.[T] We perform a detailed 
analysis of transport and radiation properties for these 
cases. 



II. GENERAL RELATIONS 

We consider a stack of intrinsic Josephson junctions 
(mesa) located at < x < L, with modulated JCC 
jj{x) = g{x)jj, where jj is the JCC density at the ref- 
erence point at which g(x) = 1. When all junctions os- 
cillate in-phase, the c-axis homogeneous phase obeys the 
following reduced equation 



ay 

5r 2 



dip d 2 ip 
v c —+g{x)smp-^ = 0, (2) 



in which the unit of length is the c-axis London penetra- 
tion depth, A c , and the unit of time is the inverse plasma 
frequency, l/u> p . Both A c and ui p are also defined at the 
reference point at which g[x) — 1. We will use these 
reduced units throughout the paper, converting to real 
units only in some important final results. The reduced 
damping parameter, v c , is related to the quasiparticle 
tunneling conductivity, v c = ATt/j c / e c oj p . We will neglect 
inhomogeneity in the dissipation parameter v c , because 
dissipation plays a minor role in the following considera- 
tion. 

In the resistive state the phase is given by 

tp = ujt + 6(t, x), 6(t, x) = Re [9 w (x) exp(— iur)] . (3) 

where w = lo /lo p is the reduced Josephson frequency. We 
will use the linear approximation for the oscillating phase 
9(t) valid for 9{t) < 1. As sin(o)r) = Re[z exp(— iQt)], 
the amplitude of the oscillating phase, 9^, obeys the fol- 
lowing equation 
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The boundary conditions follow from the relation be- 
tween the oscillating electric and magnetic fields in out- 
side medium and the Josephson relations between the os- 
cillating phase and these fields^. In general, the bound- 
ary conditions for the c-axis homogeneous oscillating 
phase describing radiation can be written as 
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dx 



/T c\n 
dr'Pir ~t') — , for x = 0, L (5) 
-oo dr 



or, in Fourier representation, 

d9 u /dx = =FiC#wj f° r x = 0, L 

/>oo 

with £ = cj(3 u and /3 W = / cLt[3{t) exp(iojr). 
Jo 
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Here the kernel (3(t — r') depends on electromagnetic 
properties of the outside media. These boundary con- 
ditions assume that there is only outcoming waves at 



both boundaries. This means that we neglect reflected 
waves, coming back to the stack and mixture of radia- 
tion coming from the opposite sides. Such a mixture can 
be suppressed, if the mesa is bounded by large metallic 
contacts on both sides acting like screens. Derivation of 
the boundary condition in such a situation is presented 
in Appendix |A| and gives the following result for £(^0 
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The radiation losses are determined by the real part of 
£. Its imaginary part only slightly displaces the reso- 
nance frequencies. In the following, we will neglect the 
imaginary part in the analytical estimates. Using typ- 
ical values N = 1000 and A c = 185/imP, we estimate 
sN/2e c X c = 3.5 • 10~ 4 , indicating that |£| is typically 
very small. 

The oscillating phase determines transport and radia- 
tion properties of the mesa. Without interference, the 
total radiation loss, Ptot, is a sum of radiation pow- 
ers coming from the left and right, sides, P tot (w) = 
fleft(^) + -Fright (w). The left-side power, Pi eft , is given by 
Poynting vector at the boundary, which is determined by 
the oscillating electric and magnetic fields at this side. 
These fields, in turn, may be related to the boundary 
value of the phase, 0^(0). In the case of the boundary 
with free space, in real units, FWt(w) can be presented 

oJUl 
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Correspondingly, the power radiated from the right side, 
-Fright) is obtained from this formula by replacement 
6U0) -► 9 U {L). For uj/2i: = 1 THz and = 1000, 
we obtain an estimate for the prefactor, 
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This estimate provides the upper limit for possible radi- 
ation power in the case of strong resonance 9 U ~ 1. 

We can also obtain a useful general expression for the 
power conversion efficiency, Q = Ptot/ (jzE z L), the frac- 
tion of the total power fed to the junction which is con- 
verted to radiation. The total power fed into the stack 
can be represented as 



J Z E Z = 



2 7r(47rA c ) : 



: {v c w + ij), 
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where v c co is the quasiparticle current a c E z in units of jj 
and ij = (g(x) sin (j>) is the reduced JCC density, 
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Combining Eqs. Q and (|9}, we derive 
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FIG. 1: (Color online) Mesas with modulation of the Josephson critical current density considered in the paper: (a) linear 
modulation (b) parabolic modulation, and (c) steplike suppression of the critical current near the edge. The lower plots illustrate 
the shapes of lowest excited Fiske resonance modes. 



In the following sections, we will first consider the be- 
havior near the resonance frequencies and then we will 
present the most interesting special cases of modulation 
allowing for exact solutions. 



boundary conditions ^ these terms can be transformed 
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ill. ENERGY BALANCE 

We consider first the energy-balance relations. The 
reduced energy in units of L z L y ^/ (2>2tt z s 2 \ c ) accumu- 
lated in the phase oscillations inside the stack is given 

by 



dx 



1 see 

2 I dr 
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Therefore, the energy-change rate is given by 
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As follows from the Eq. ([2]), the first term can be trans- 
formed to 



dx- 



86 


'8 2 9 


8 2 9 


dr 




dx 2 




dx ^ 


89\ 2 


f 




Jo 


8t) 



/ dx— g[x) sinwr. 
Jo ot 



Here, the first term accounts for the quasiparticlc damp- 
ing, while the second term gives the driving force from 
the Josephson oscillations leading to pumping of energy 
from a dc source into the electromagnetic oscillations in- 
side the stack. The last two terms in Eq. ( 13 I account for 
the radiation losses at the boundaries. For the general 



where we introduce a notation for the operator /3^p = 
dr' P(t — t') Jp. Therefore, the total energy-change 
rate can be written as 
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For a steady state, the energy has to remain constant, 
meaning that the energy supplied by the Josephson os- 
cillations has to be exactly compensated by the quasipar- 
ticlc and radiation losses. 



IV. BEHAVIOR NEAR RESONANCES FOR 
ARBITRARY MODULATION 



In this section, we obtain approximate results for the 
current and radiation in the vicinity of the resonance 
frequency u) m = mm /L, where the phase can be approxi- 
mated as the corresponding cavity mode 



6{x) w ^cos(rmrx/L). 



(15) 



We neglect the small influence of the radiation on the 
shape of the resonance mode. In this case, the energy 
in the mode (12) and energy-change rate (14) can be 
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approximated as 
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is the coupling parameter. Therefore, equation for the 
mode amplitude is given by 



d 2 ip _ 2 , dip 
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Using complex representation ip = He[ipuj exp(— Hot)], we 
obtain a solution 



&m + i i V c + V r ) Co ' 



where 
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L 
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p 



is the parameter of the radiation damping (the last for- 
mula is written in real units). One can see that both the 
quasiparticle dissipation and radiation contribute to the 
resonance damping P3 The cavity quality factor is given 
by Qc — u m j(y c + v r ). Optimal power conversion is 
achieved when the main contribution to damping is com- 
ing from the radiation, v c <C v r . Comparing the damping 
channels using Eq. ([7]), we obtain that this is achieved 
for a sufficiently large number of junctions in the stack 



N>N n 



e c v c L 2ira c L 
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Taking s = 1.56nm, we also rewrite this formula in the 
practically convenient form as 



N„ 



576 o- c [l/(fi • cm)]L[/H//[THz] 



For typical values a c = 0.003-0.01 (fi-cm)" 1 , L ~ 40A«n, 
and / = to/2ir = 1 THz, we estimate N a = 70 - 250. 
In the regime of dominating radiation losses, the cavity 
quality factor is simply given by Q c = e c L/(2L z ). 

The solution (20 1 allows us to obtain the average JCC 



ij = (g(x) sin {lot + Re \tp w exp(— Hot)] cos(m7ra/L)}) 

1 dli i V c + V r ) LO 
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f + {v c + Vrf Co 2 



(23) 



This gives the maximum current enhancement in the res- 
onance 

3^/4 



ij.i 



(v c + V r ) Ur, 



(24) 



A similar result has been derived in Ref. |2H for the case 
of a single junction without radiation losses. Comparing 
this result with the reduced quasiparticle current, v c ui, 
we see that the resona nce feature i n I-V characteristic is 
pronounced if g m > 2\J (y c + u r ) u c Co. In the case of dom- 
inating radiation losses we can rewrite this condition in a 
more transparent form, g m > 2yj2L z v c j '(e c L)Co 3 / 2 . For 
v c = 0.01, Co = 10, L z = 1.5/im, and L = 40/zm corre- 
sponding to m = 1, we obtain that the resonance feature 
in the I-V dependence becomes strong when the coupling 
parameter exceeds 0.5. In the case of strong resonance, 
the total current i(Co) — u c u> + ij nonmonotonically de- 
pends on the Josephson frequency Co (voltage). In this 
case, only the increasing part di/dCo > is stable. 
The total radiated power from both sides is given by 

P tot (co) = 2P SC ^(3^ = 2P sc Re[f3 UJ ]Cb 2 \iP UJ \ 2 



2P sc f3 u] w 2 g 2 



\t , j 2 — /, )2 



+ (v c + v r ) Co 2 



(25) 



Here, the scale of P tot is given by P sc = L y L z X c E p jj/2, 
where E p = ^olo p /(2ttcs) is the electric field correspond- 
ing to the plasma frequency. For the maximum ra- 
diated power in the resonance we obtain Ptot(w m ) = 
2P sc f3 UJ g 2 n / (y c + v r ) . In the regime of dominating ra- 
diation losses, v c <C v r , using Eq. ^ and jj = 
c$>o/(8ir 2 s\ 2 ), we obtain a very simple and universal es- 
timate for maximum total radiated power (in real units) 



Poti^m) 



2lo 



(26) 



An important observation is that for a tall stack in res- 
onance, the radiated power does not depend on N due to 
compensation between the factor A 2 in front of the ra- 
diated power ^ and the amplitude of phase oscillations 
in the resonance, which, due to the increasing radiation 
losses, drops as £ -2 cx A -2 , see also Ref. [TjJj This com- 
pensation only exists in the regime when the damping of 
the resonance is caused by the radiation which is realized 
under the condition ([22 1 . 



The power conversion efficiency is given by 



Q 



P, 



L (ij + v c Co) Co 



(27) 



In resonance, it can be represented in a quite transparent 
form as a product of two factors 



g 2 n + 4 (u c + v r ) v c Cb 2 v c + v r 



(28) 



where the first factor represents the relative current in- 
crease in the resonance, ij ima ,x/(v c Co + «j, m ax), and the 
second factor is the relative contribution of the radiation 
to the resonance damping. We can see that, remarkably, 
in resonance the conversion efficiency can approach 100% 
provided (i) the resonance feature is pronounced in I-V 
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Long mesa, kj_ y > 
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V r ~ 2v r o 


V r « Vh)[1 + •/<)(&„,£)] 


v r « 2v r0 [l + Jo (WO] 



FIG. 2: (Color online) The parameter of radiation damping for the fundamental mode for different mesa designs in the case 
k^L z <C 1: long mesa with screens mostly considered in the paper (left), rectangular capacitor (middle), and mesa with ground 
plate (e.g., mesa fabr icat ed on the top of bulk crystal) (right). Here v r o = uL z / (s c lj p L) and the function Ti t o(a m , a y ) is defined 
in Appendix IB] Eq. (B8l. In the regime of dominating radiation losses, the cavity quality factor Q c is directly determined by 

V r aS Q c = UjJ(l/ r OJp) . 



dependence, g m > 2yJ(v c + v r ) v c lo and (ii) the losses are 
dominated by the radiation, v c < v r . Both conditions are 
quite realistic. 

Remember that simple and transparent results for the 



typical number of junctions (22 1, current (24 1, radiation 



power ( 26 1 , and conversion efficiency ( 28 1 are valid only 



in the case ujL v /c > 1 and without mixing of radia- 
tion coming from the opposite sides. These results can 
be generalized for other cases. Radiation losses of the 
resonance mode in the short rectangular mesa can be 
approximately calculated similarly to radiation out of a 
rectangular capacitor.^ These calculations are summa- 
rized in Appendix [S] and the results for the radiation 
damping parameter for different cases are presented in 
Fig. [2] In the case of long mesas, k w L y ^> 1, the radi- 
ation damping parameters for different geometries differ 
only by numerical factors of order unity. In the opposite 
limit, k u L y -C 1, v r acquires an additional small factor 



Linear modulation 



In this section, we consider a mesa, with linearly mod- 
ulated JCC, g(x) = 1 — 2rx/L, see Fig. [T^l. Here, we 
introduce the new coordinate, x = x — L/2, which is 
symmetrical with respect to the mesa, —L/2 < x < L/2. 
This means that the JCC at the left side is larger by the 
factor (1 + r)/(l — r) than the JCC at the right side. 
Such a modulation couples homogeneous Josephson os- 
cillations with the odd Fiske modes Q, m = 21 + 1, 
including the fundamental mode, m — 1, and the cou- 
pling parameter (18 1 to this mode is connected with the 



modulation parameter as 



.91 = 8r/w 2 



(29) 



SPECIAL CASES OF MODULATION 



In this section we consider several special cases of mod- 
ulation for which the problem allows for exact analytical 
solution and full analysis of transport and radiation prop- 
erties. We will consider three cases: linear modulation, 
parabolic modulations, and steplike suppression of the 
critical current. Practical ways to prepare such modula- 
tions are suggested in the discussion section ??. 



Solution of Eq. Q with the boundary conditions ^ in 
the case of linear modulation can be found exactly. Split- 
ting the solution into the symmetric and antisymmetric 



7 




FIG. 3: (Color online) Representative plots of the Josephson-frequency (or voltage) dependences of (a) the current density j, (b) 
radiated power from the left side, P, (c) power conversion efficiency Q, and (d) amplitude of oscillating phase at the boundary 
for stacks with linearly modulated JCC with parameters r — 0.2 and r = 0.4. The unit of current densi ty i s the JCC density, 
jj, the frequency unit is the plasma frequency u> p , and the unit of radiated power is 10 _3 Pj, see Eq. (36 1, corresponding to 
P/Ly ~ 0.16W/cm. The following parameters have been used L = 0.23, u c = 0.006, and Re[£] = 0.00035£ 2 . For comparison, 
the case of homogeneous JCC (r = 0) is also shown but P, Q, and \9\ 2 are indistinguishable from zero at the scales of the plots. 
In this case, P « 1.1 • 10~ 4 and Q « 10" 4 at lu = 13.4. 




FIG. 4: (Color online) Comparison between the powers radi- 
ated to the two opposite sides of the mesa for r = 0.2 and 
the same parameters as in the previous figure. One can see 
that the difference is rather small and amounts to a slightly 
different asymmetry of the peaks. 



we derive 



lirxjL 
up + iv c oj 



(2i/L + Qrsia{p u x) 



(a) 2 + iv c u) [p w cos(x) - <sin(x)] ' 
i 



Hi 1 + ii/ c u> 



(cosip^x) 



(a) 2 + iv c u)) (p u sin(x) + <cos(x)) 
with p 1 ^ = £j 2 + w c lu and x = PuL/2. 



(30) 



(31) 
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FIG. 5: (Color online) Evolution of the radiation power and 
resonance feature in the I-V dependence with increasing num- 
ber of junctions in the stack, N, for r = 0.2 and the same 
parameters as in Fig. [3] Note that while the radiation power 
increases with N, the resonant feature in current becomes less 
pronounced due to the increasing radiation damping oc iV 2 . 
The inset in the upper plot shows the iV dependence of the 
maximum radiation power, and the dashed line in this plot 



indicates the large- iV limit obtained from Eq. ( 26 1 
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Only the antisymmetric phase is coupled to the resonance 
mode. In particular, for the boundary phases, we have 



L 

±- 
2 



± 



±- 



L 



ir [-cos(x) +sin(x)/x] 
Pu> [Puj cos(x)-iCsin(x)]' 

^sin(x) 

[p w sin(x)+<cos(x)]' 



(32) 
(33) 



The radiated power is determined by the boundary 
phases by Eq. ([8]), where we have to replace 9 W (0) with 
6l a) (±L/2)+e¥(±L/2). 

Near the resonance ujL — it, using p^ w u) + iv c /2 and 
cos(x) ~ (tt/2 — QL/2) — iv c L/A, we obtain 



8ir/Tr 2 



(u + iv c ) [2(a) — 7r/Z) + i(z/ c + v r )\ ' 



0£°>(L/2) 



where v r is defined in Eq. pi) . Using the coupling pa- 
rameter ( 29 1 , we can see that this result is consistent 



with the general formula (20 1 near the resonance. The 



maximum antisymmetric phase in the resonance can be 
estimated as 9\y{L/2) w — 2r / (it (i^uj 2 ) . It exceeds the 
nonresonant symmetric part approximately by the factor 
r/flu ~ re c L/L z . In resonance, using Eqs. (26) and (29 1, 



we obtain for the radiation power from one side in real 
units 



16r 2 L 2 j 2 



It exceeds nonresonant emission by the factor 
(re c L/L z ) 2 . In particular, for uj/2-k = 1 THz, jj = 500 
A/cm 2 , and r = 0.4, we estimate P r /L y ps 0.05 W/cm. 

In the average JCC (10 1, the contributions from sym- 
metric and antisymmetric parts split, ij = ij. a + ij, s - 
Direct calculation gives 



iJ,a- 2 



+ Re 



v c /3 



(bJ 2 + V 2 ) UJ 

[cos(x) - sin(x)/x] (2i/i + C) Ix 



:.2 



> j.. 



= -Re 



ivJZ)) \p u cos(x) - <sin(x)] . 
i 2i^sin(x)/i 



1- 



+ iv c G)\ p w sin(x)+iCcos(x) 



(34) 



(35) 



From the general formula for ij a near the resonance we 
obtain a much simpler result 



16r 2 (z/ c + i/ r ) /u) 

^ A{Cj-^/L) 2 + {y c + u r f 



and the maximum current enhancement in the resonance 
is given by 



lj,a 



16r 2 /7r 4 
uj [y c + v r ) 



These results are also consistent with the corresponding 
general formulas ( 23 1 and (24 1 if we use the coupling 



parameter (29 1. For comparison, the symmetric part of 



the JCC at the resonance frequency can be estimated as 



',/... 



2w 3 



As expected, the resonant enhancement of the current 
exceeds the nonresonant radiation correction by the same 
factor (r//3 w ) 2 = (re c L/L z ) 2 as for the radiation power. 

For illustration, we present the behavior near the reso- 
nance for mesas with two modulation parameters, r = 0.2 
and r = 0.4. As a unit of the radiation power, we selected 
the quantity 



Pj = LyXJj/UJp, 



(36) 



which is independent on the sizes L and L z . This choice 
of unit is suggested by the result p6] ). For A c = 185 
/xm Pj/Ly fa 163 W/cm. Figure [3] shows the Josephson- 
frequency dependences of (i) the current density j (in 
units of the JCC density in the center), (ii) radiated 
power P (in units of 10~ 3 Pj), (iii) the power conver- 
sion efficiency, Q, and (iv) the amplitude of oscillating 
phase at the boundary. For comparison, the case of ho- 
mogeneous mesa (r = 0) is also shown. In calculation we 
used the following parameters: L — 0.23, v c — 0.006, 
and Re[C] = 0.00035w 2 (corresponding to N « 1000, 
A c w 185 fim, and <r c « 0.003 [ft cm] -1 ). We can see 
that the modulation leads to the appearance of a strong 
resonance feature in the I-V dependence. Note that only 
I-V regions with positive differential resistivity are stable. 
Current enhancement in the resonance is mainly caused 
by the generation of the powerful electromagnetic wave 
and it is accompanied by a huge enhancement of outside 
radiation. The maximum radiation power for used pa- 
rameters for the case r — 0.4 corresponds to ~ 0.05W/cm 
and it exceeds the nonresonant radiation from the homo- 
geneous mesa by more than 3 orders of magnitude. It is 
important to note that the power conversion efficiency is 
also strongly enhanced in the resonance, reaching 20% 
for r — 0.4. The plot of \8\ 2 shows that for selected pa- 
rameters it remains smaller than one in the resonance 
and, therefore, the linear approximation used in calcula- 
tions is not violated. 

In spite of the asymmetry of the JCC, the powers ra- 
diated to the opposite sides of the mesa in resonance are 
approximately the same, because the radiation is mostly 
promoted by the resonance mode which has identical am- 
plitudes of the oscillating electric field at the opposite 
sides. This is illustrated in Fig. [4] where these powers 
are plotted for r = 0.2. One can see that the peaks have 
slightly different asymmetries originating from the sym- 
metric phase. 

Figure[5]illustrates the evolution of the radiation power 
and resonance feature in the I-V dependence with in- 
creasing number of junctions in the stack, N, for r = 0.2 
and the same parameters as in Fig. [3] The number of 
junctions above which the radiation losses dominate (22 1 
can be estimated for used parameters as N a w 75. We 
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can see that the current and radiation have opposite ten- 
dencies: while the radiation power increases with N, the 
resonant feature in current becomes less pronounced due 
to the increasing radiation damping. 



B. Symmetric parabolic modulation 

In this section we consider a symmetric modulation. 
For simplicity, we assume a simple parabolic profile of 
the JCC density g(x) — 1— r(2x/L) 2 , see Fig.[l3, where, 
again, x — x — L/2 and 1 — r is the ratio of JCCs at the 
edge and in the center. The cases r > and r < 
correspond to current suppression and enhancement at 
the edges respectively. Such a modulation will lead to 
excitation of only even frequency modes ([lj, m = 21. 
In the following, we will focus on the lowest even mode 
with m = 2. To have this resonance at LJ2/2ir = 1 THz 
assuming e c w 12, the mesa size has to be rather large, 
L = 86.5 (im. The coupling parameter p8| to this mode 
in our case is given by 



.92 = -ir/n 2 



The oscillating phase in the case of parabolic modu- 
lation also can be found exactly. From symmetry, the 
solution of Eq. Q with boundary conditions ^ must be 
an even function of x, and it has the following form 



2ir(2/L) 2 
(u 2 + iv c iv) 2 



i(l-r(2x/L) 2 ) „ , . /nn , 
- -o , - +Ccos( Pul x), (38) 



where the first two terms give a particular solution of the 
inhomogeneous equation and the last term is the solution 
of the homogeneous equations. From the boundary con- 
ditions ([6|, we obtain 



C 



-Ur/L + C (2r/x 2 + 1 - r) 
[Pu sin(x) + iCcos(x)] (w 2 + iv c uS) ' 



(39) 



where and x are defined after Eq. (|31| . The oscillating 
phase given by Eqs. (38) and (39 1 determines all other 
observable properties. 

The powers radiated from both sides (TSJ) are deter- 
mined by the boundary phase W) o = Oj{h — ±L/2), 



which we obtain from Eqs. (|38| and (|39|), 

2ir(2/L) 2 



i(l - r) 



(lj 2 + ivjb) u + 
[-ter/L + C (2r/x 2 + 1 - r)] cos(x) 
{Pu> sin(x) + K cos(x)) (w 2 + ivjb) 

In the resonance, uiL = 2n, using 9 u (L/2) » — 2r/(n(uj), 
we obtain for the power radiated from one side in real 
units 



±L y L 2 r 2 j 2 



(40) 



This result is consistent with the general formula (26 1 
with the coupling parameter (37 1. For u/2ir 
j c = 500A/cm 2 and r 



1 THz, 
we obtain an esti- 



= ±0.5, 

mate for the radiated power in the resonance, P/L y f» 
0.085W/cm. Note that jj in Eq. (go) is the JCC density 
in the center while the the JCC density at the edge, jj >e , 
is given by jj e = (1 —r)jj. In the case of r < vanishing 
of superconductivity in the middle, which corresponds to 
the limits jj — > and r — > — oo, does not lead to van- 
ishing of radiation power because r 2 j 2 — > j 2 e and the 
radiation is determined by the JCC density at the edge. 

To find the voltage-current characteristic, we calculate 
the average reduced JCC (10). Using oscillating phase 
(38 1, we obtain 



2Cb (Q 2 
C 
~2 



Re 



sinx 
X 



2r 

y 



-2r 



r 2 


16r/L 2 


y + 


UJ 2 + V 2 


2 N 


\ sin X 




' 2x + 



cosx 



(41) 



(37) At the resonance frequency, ujL = 27r, we estimate 



e c r' 



! L 3 /(27r 6 L z ) 



Figure [6] shows the representative Josephson-frequency 
dependences of the current density j, radiated power P 
to one side, the power conversion efficiency, Q 7 and the 
amplitude of oscillating phase at the boundary for neg- 
ative modulation parameters, r = —0.2 and r = —0.4, 
corresponding to the case of stronger superconductivity 
at the edges. In calculation, we used the following pa- 
rameters: L = 0.8, v c = 0.005, and Re[C] = 0.003a; 2 
(corresponding to N w 1000 and A c w 200/im in the cen- 
ter). Overall, the behavior is very similar to the case of 
linear modulation shown in Fig. [3] with minor quantita- 
tive differences. We also see that the modulation leads to 
the appearance of a strong resonance feature in the I-V 
dependence accompanied by a huge enhancement of the 
outside radiation and power conversion efficiency. 



C. Steplike suppression of critical current near the 

edge 

In this section, we consider the case when there is a 
region with suppressed JCC on one side, see Fig. [lj; 

f 1, for <x < L- W 
^ = {l-r, i orL-W< X <L- ^ 



The coupling parameter ( 18 1 to the fundamental mode 
in this case is given by 



2r 



fji 



sin(7rW/L). 



(43) 



The solution of equation for the oscillating phase Q 
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FIG. 6: (Color online) The Josephson-frequency (or voltage) dependences of (a) the current density j, (b) radiated power 
to one side P, (c) power conversion efficiency Q, and (d) amplitude of the oscillating phase at the boundary for stacks with 
parabolic profiles of the JCC and negative modulation parameters, r = —0.2 and r = —0.4, corresponding to the case of current 
enhancement at the edges. All units are the same as in Fig. [3] The following parameters have been used L = 0.8, v c — 0.005, 
and Re[£] = 0.003iD 2 , corresponding to N w 1000 and A c w 200/im in the middle. For comparison, the case of a homogeneous 
mesa (r = 0) is also shown. In this case, P w 5 • 10~ 5 and Q w 5 ■ 10~ 5 at uj = 7.85 
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FIG. 7: (Color online) The radiated power and I-V depen- 
dence near the resonance for different widths of the suppressed 
region near the mesa right side. The JCC density in the sup- 
pressed region is assumed to be half of its value in the rest 
part, r — 0.5. All units and parameters are the same as in 
Fig-i 



can be built in the piecewise form, 



i/(uj 2 +iv c Cj) + A + exp(zp w x) 
+A- exp(— ip u x), for < x < L — W 

i(l - r)/{tb 2 + iu c u)) + (A + — C+) exp^p^x) 
+ (A_-C_)exp(-ip w a;), for L-W <x<L 



(44) 



with plj = uj 2 + iv c uj. Matching U and d9 u /dx at x 
L — W , we obtain 



-ir/2 



uj z + ii> c uj 



- exp [Tipu, (L - W)} 



(45) 



Using this result, from the boundary conditions ^ we 
find the coefficients A± , 
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A _ ({(Pu, ± Qexp(=Rx) + {Pu T C) [1 ~ r(l - cost?)]} - ir PuJ ( PuJ =p C) sinr? 
2 (lj 2 + iv c Lj){[pl + C 2 } sin x + Zip^C cos x} 

where x = an d f] = Pu>W . This gives for the boundary phases which determine the outside radiation 

* , C{Pu [1 + cosx] - iCsinx} - rpu{£(l - cosr?) +ip ul siarj} 



^(0) : 
B U {L): 



Ul^ + lV r LO 



(ui 2 + iv c Cb) [(pi + C 2 ) sin x + Zip^C cos x] 



i {1 — r(l — cos 77)} 
a) 2 + iz/ c w 

C [pu, ( 1 + cos x) - < s m X] ~ r [C ( 1 - cos t?) + ^ sin (??)] [p^ cos x ~ i_C sin x] 
(w 2 + w c a>) { [pi + C 2 } sin x + 2ip UJ ( > cos x} 

I 



(46) 



(47) 



(48) 



Near the resonance, QL = ir, the coefficients A + and A_ 
can be strongly simplified 



A+ pa A_ 



-irsin[wVy]/(2w 2 



2(cj — 7r/L) — i(i/ c + f r ) 
In the resonance A + w A_ w rsin [ojW] /(4d) 2 /3 u ) giving 

i r sin [tliW] 



uj 2 + iv c Cj 2uj(3^ 



At W <C i the resonance is strong, W (O) > 1, if rW ^> 
2d)/3 tlJ or, in real units, 



rW > 



t t 2 A 2 L, 
e c i 2 



as 



which is consistent with the general formula ( 26 1 



The reduced JCC flowing through the stack ( 10 1 can 
be computed as 



. _ {L-r{2-r)W)v c (l-r)r ^ 



2L (u 2 + v 2 )Q 
Re 



2L 

M/x 



sin 7? 



(Q 2 + iv c u) { [pi + C 2 ] sin x+ 2ip L0 C cos x} 



(50) 



The total radiated power in resonance can be estimated with 



[Puj COS - 
, ■ 2 • 2 V 

+ ir sin — 



i£sin ■ 



{Csh 



V ., ■ V 

cos iC sin — 

2 s 2 



p w cos 



Ccos- 

. _ V 
2 



I V V\ XI 1 

+ i [p u cos - - < sin - J sin - j 



+ i( sin 



Near the resonance, we estimate JV w 
—ip 2 r 2 sin 2 (tDW 7 ) /2 which gives 

. r 2 sin 2 (^)K£/2 + 2^}/2tt 
w 2 |sin 2 (£L) + (y c Lj2 + 2/3 w ) 2 | 

The maximum current enhancement in the resonance is 
given by 



- 2 sin 2 (irW/L) 
7r 2 o> (u c + ApJL) ' 



General cumbersome formulas (47 1, (48 1, and (50 1 can 



be significantly simplified if we assume the conditions v ci 
ICI <C uj and W <C L valid in most practical situations. 
In this case these equations can be represented in approx- 
imate, simpler form, 



/q\ _ i irW/ti 

' ub 2 + iv c uj sin x — i (v c L/2+2{3 UJ ) cosx' 



irW cos x/w 



cj 2 + w c tl) sinx — i {v c L/2 + 2(j u ) cosx 



, (53) 



in agreement with Eqs. (|24|) and (43 1 



and 
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v c v c rW 

lj fts V c LO 



Re 



2u 3 2Lti 3 

j 2/3 u sinx - iurW [sinx - 2i/3 OJ [cos (x - 5) + cos (§)] - 2wrWcos (x - f)] 
I 2L<I> 3 [sin x - « (v c L/2 + 2/3 u ) cos x] 



with x ~ and 77 « £>W\ 

Figure [7] illustrates the evolution of the radiated power 
and resonance feature in the I-V dependence with in- 
creasing width of the suppressed region. The JCC den- 
sity in the suppressed region is assumed to be half of 
its value in the rest part, r — 0.5. For used parame- 
ters, the maximum radiation power in this plot is around, 
Pjxux/Ly ~ 0.05W/cm. 

VI. DISCUSSION AND SUMMARY 

Let us discuss now practical ways to prepare mesas 
with lateral modulation of the critical current density. 
Mesas with linear JCC modulation can be fabricated in 
a crystal with inhomogeneous doping. One possible way 
to prepare such inhomogeneity is to utilize the sensitiv- 
ity of doping in BSCCO to the oxygen concentration. 
Due to strong temperature dependence of the oxygen 
diffusivitjESI, in principle, the oxygen concentration pro- 
file in the crystal can be prepared by short-time anneal- 
ing by carefully selecting the annealing temperature and 
time. In a similar way, mesas with parabolic-like profiles 
can be prepared by short-time annealing of mesas them- 
selves already after fabrication. Another way to prepare 
modulation in a controlled way is to use radiation with 
high-energy electrons, protons, or heavy ions. If part 
of the mesa is protected by a mask, this radiation will 
produce a mesa with steplike suppression of the critical 
current. 

The major technical challenge is to prepare a mesa 
with significant modulation of the Josephson coupling 
identical in all junctions. Variation of parameters in dif- 
ferent junctions, which may be caused by composition 
variations, inhomogeneous heating, and different junc- 
tion areas would strongly reduce the optimal perfor- 
mance. The quantitative analysis of the radiation prop- 
erties of mesas with such parameter variations in different 
layers will be done elsewhere. 

As the optimal mesa size is rather large, another major 
technical problem is sample heating due to quasiparti- 
cle damping. The self-heating in the BSCCO mesas has 
been investigated by many experimental groups^. The 
major focus of these studies was the influence of heat- 
ing on the gap feature in I-V characteristics, which is 
located at voltages 30-60 mV/junction. Even though our 
voltage range is significantly lower, ^2 mV/junction, the 
heating is still expected to be significant due to the re- 
quired large lateral size of the mesa. For example, for 



a c = 0.003 cm], N — 1000, and L y = 300Atm, 10 

mW of power will be dissipated inside the mesa. This 
heat has to be removed from the mesa faces. Therefore, 
efficient heat removal is crucial for operation of the de- 
vice. Recent experimental observations of the resonant 
emission using underdoped BSCCO^Sl demonstrate that 
the heating effects can be manageable even in large-size 
mesas with lateral sizes of several hundred micrometers 
in the voltage range corresponding to the Josephson fre- 
quencies around 1 THz. 

The designs with improved thermal management may 
include, for example, fabrication of underdoped mesas on 
the top of overdoped crystal, using massive gold contacts 
on the top and bottom of the mesa, and placing an insu- 
lator with high thermal conductivity, such as sapphire, in 
good thermal contact at the side of the mesa. From these 
considerations, mesas with asymmetric modulation look 
more preferable than ones with symmetric modulation, 
because they need a smaller lateral size for the lowest 
resonance mode. In the case of symmetric modulation, 
the design with suppression of the JCC in the middle, 
r < 0, looks more practical for better thermal manage- 
ment. In fact, the material in the middle can even be 
made insulating, because this part is needed only to form 
almost standing wave at the working frequency. To excite 
resonance mode, it is sufficient to have superconducting 
regions only at the edges. 

In conclusion, we demonstrated that a stack of the in- 
trinsic Josephson junctions with modulated Josephson 
coupling represents a very powerful and efficient source of 
electromagnetic radiation at the resonance frequency set 
by its lateral size. Selecting this size, the generation fre- 
quency can be tuned to the terahertz range. Power levels 
up to several milliwatts look plausible in such structures. 
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APPENDIX A: BOUNDARY CONDITIONS FOR 
THE HOMOGENEOUS OSCILLATING PHASE 

In this appendix, we consider the boundary conditions 
for the oscillating phase at the edges and the radiation 
power for a stack of intrinsic Josephson junctions. We 
will limit ourselves to the case when the oscillating phase 
is identical in all junctions. A more general case will be 
considered elsewhere. The oscillating phase 9^ defined 
by Eq. ^ is connected with the electric and magnetic 
fields by the Josephson relations 



-I 

2ttcs 
$ 



(Al) 
(A2) 



Therefore, the boundary conditions for the oscillating 
phase at the edges are determined by the relation be- 
tween the fields E z and B y in the outside media, which 
we assume to be monochromatic with time dependences 
oc exp(-iujt) . 

Outside dielectric media at \x — L/2\ > L/2 is charac- 
terized by the dielectric constant e d , and we assume only 
outgoing wave in this space. The Fourier components of 
fields with \k z \ < y/e^k^ propagate in the media, while 
the field components with \k z \ > y/e^k^ decay. In par- 
ticular, for E z (u),x,k z ) at x > L, we have 

E z (uj,x,k z ) = E z (u>, L, k z ) exp [ik x (u, k z )(x - L)] (A3) 

with 



yje d kl - fc2sign(w), for \k z \ < y/ed\k u \ 



isjk 2 z - e d kl, for \k z \>^/s d ~\k 



k x (uj, k z ) - 

I i^K%-e d KZ, ior \k z \ > y/e d \K U \. 

(A4) 

Other field components, E x and B y , are also expressed 
via E z (lo^ L, k z ). First, E x (oj, x, k z ) is obtained from Eq. 
(A3) and the Maxwell equation V • E = 0, and then 
B v (lu, x, k z ) is obtained using the Maxwell equation (V x 
E) y = ik^By leading to the following result 

B y (uj,x,k z ) = -E z (uj,L,k z )— — ^— exp [ikAw,k z ) (x-L)] . 

(A5) 

This gives the relation between the fields at the boundary 

x = L, 



By(L,k z )-- 

(((J,k z ) = 



-({u,k z )E z (L,k 
\e d /y/s d kl - 



(A6) 



-ik UJ e d /^/W z 



- fcf, for \k z \ < y/e^lk^ 
e d k1, for \k z \ > Je2\K 



The condition at x — has opposite sign, B y {Q,k z ) — 
((uj,k z )E z (0,k z ). Note again that the term ((oj, k z ) 
for \k z \ < y/Sdk u originates from outgoing electro- 
magnetic wave (radiation), while the term £(w, fc z ) for 
\k z \ > y/^dk^j is due to the wave decaying at distance 



term does not carry energy out of the junctions. In par- 
ticular, for k z — 0, we have B y (L,0) — — y/e^E^L^ 0), 
leading to the simple boundary condition for the ho- 
mogeneous oscillating phase in the limit L yi L z 3> Aq, 



Va^o, = ^(i^fesxi I c)Q^ for x = L, 0. The relation (A6l 



can also be rewritten in the frequency-space representa- 
tion as 

/oo 
dz'U(z- z',w)E z (L,z' i uj), (A7) 
- OO 

U(z,u) = -y [IfcJJotv^l^D + ik u N (^/e^\k u z\)] , 

where Jq(z) and No(z) are the Bessel functions. 

The same approach can be used in the realistic case 
of a crystal small along the z axis, L z < X Q , if we know 
the radiated electric field at the planes x — 0, L outside 
of the crystal, at \z\ > L z /2. If we put well conducting 
screens there, we can approximate E z — at \z\ > L z /2. 
In this case for the homogeneous n-independent electric 
field, we obtain for the average magnetic field at the edge 



L z e d 



2i 

\kj\ feu hi 



C 



E z (L,cj), 



7T \f£d,L z \k u 

(A8) 

with C = 2exp(3/2 - lE ) w 5.03, where j E w 0.5772 is 
the Euler constant. We can sec that for a small-size mesa, 
the magnetic field at the boundary is reduced by the 
factor ~ L z k u in comparison to the infinite-L z case. This 
gives the following boundary condition for the oscillating 
phase 



ik u L z e d 



kuj k^j In — 

7T y/e d L z \k u 



(A9) 



kt) 1//2 from the crystal boundary. The latter 



for x = L,0. This corresponds to the boundary condi- 
tions (JsJ) and Q in reduced coordinates and e d = 1 used 
in the paper. Therefore, for short crystals, the boundary 
condition can not be written in the form of an instan- 
taneous relation in between the space and time deriva- 
tives of the phase. This significantly complicates their 
numerical implementation. Screens also completely iso- 
late semispaces x > L and x < and eliminate interfer- 
ence of radiation coming from the opposite edges. 



APPENDIX B: RADIATED POWER FROM A 
RECTANGULAR MESA 

The radiation power from a short rectangular mesa 
can be found approximately. For such a mesa radia- 
tion influences weakly the shape of the resonance mode. 
In such a situation, the radiation is mostly determined 
by the distribution of the oscillating electric field at the 
mesa edge, which, in turn, is determined by the shape of 
the internal mode. Finding the radiation occurs to be a 
somewhat easier problem than finding general boundary 
conditions for the oscillating phase. An approximate ex- 
pression for the radiated electric field far away from the 
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related to the oscillating phases 9 U at these edges accord- 
ing to Eq. (Al). We consider modes inside the crystal 



which arc uniform along the z-axis (synchronized Joseph- 
son oscillations in all IJJs). In this approximation, the 
electric field E z inside the crystal is described by standing 
waves characterized by the indices m and n: 

E z (r) = E z (m, n) cos(k x , m x) cos(k y ^ n y), (B2) 

k Xy m = TfrTl/ 'L Xl ky n = TTTlfL/y 

Faraway radiated electric field in terms of E z (n,m) is 
given by the expression 

E ^_ k u L z exp(ik„r) ! Ms(r ' )exp( _j fc ^ er)(erXer)d ^ 



FIG. 8: (Color online) Geometry for radiation out of a rect- 
angular mesa. 



crystal and, thus, the radiation power can be calculated 
using the Huyghcns' principle, as it is developed in the 
theory of antennas, see, e.g., Ref. [3T] This approach 
has been applied to resonance modes in rectangular ca- 
pacitors in Ref. HH1 Such a consideration clearly shows 
the role of screens and crystal geometry in the formation 
of the radiation. The Huyghens' principle in the formu- 
lation of Schclkunoff (equivalence principle) states that 
we can find fields outside of real sources (currents and 
charges) if we know equivalent sources placed on some 
boundary surface surrounding real sources. In particu- 
lar, the equivalent magnetic current, M s , is related to 
the tangential components of the electric field E on the 
surface as 



M s = — n x E, 
4tt 



(Bl) 



where n is the normal vector on the boundary surface. 

In the following, we consider the crystal inside the vol- 
ume < x < L x , < y < L v , and < z < L z bounded 
by thin metallic contacts on the top and bottom, see Fig. 
[8] The contacts are highly conductive, and there we can 
neglect the tangential electric field. At the crystal edges 
x = 0, L x we can neglect the magnetic field when the 
radiation power is small, as in the case of radiation from 
a capacitor with a small distance L z between electrodes, 
t u L z <C 1, see Ref. [25] In this case, we need to account 
only for the electric field at the crystal edges, which pro- 
duces the magnetic equivalent currents ( B 1 1 . They are 



where the integral is taken over the perimeter of the crys- 
tal, and the coordinate system as well as definitions of the 
unit vectors e r and By are given in Fig. |8l Integration 
over contour £ gives the following resuHpST 



E 



E z (m,ri) 



L z exp(ifc^r) 
Airr 



k h 



k 2 
x.m 



h 2 



G X Gy , 



where 



sin <p eg + cos u cos < 







COS 


• = sm6< 






_ sin 4> J 



and 



The interference factors G x = 1 — (—1)™ exp (ik^L x ) and 
G y = 1 — (— l) m exp (ik^Ly) describe the contribution of 
waves coming to the faraway point r from opposite sides 
of the crystal x' = and x' — L x along the x axis as well 
as from opposite sides y' = and y' — L y , respectively. 
Their role in the formation of the radiation becomes clear 
if we will consider different modes. From the radiated 
electric field, we can compute the total radiated power 
as 



P : 



sin^[|^ e | 2 +|^| 2 ] . (B3) 



For homogeneous oscillations m — n — 0, we obtain 



E 



3,(0,0) 



corresponding to 



E 



E z (0,0) 



L z exp(ifc^r) / k, 



Aur 



kf 



- y t [1 - exp (-ik(-L x )] [1 - exp (-ik v L y )} 

Kyi 



L z exp(ikur) [1 — exp (—ik^Lx sin 9 cos 0)] [1 — exp {—ik u L y sin 9 sin 0)] 



47rr 



sin v cos m sin < 



(B4) 
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and = 0. The radiated power (B3) is given by 
c\E z (0,0)\ 2 



P = 



32tt 3 



dO 



[1 — cos {k u L x sin 9 cos 4>)] [1 — cos {k u L y sin 9 sin </>)] 



sin f cos^ 



>sm 



For small-size crystal L x ,L y <C fc w 1 with almost uniform 
Josephson oscillations, we obtain 



D _ c\E z (0,0)\ 2 ]At2t2t2 



(B5) 



This is the result for dipole radiation because all sizes of 
the crystal are small in comparison with the wavelength 
of the radiated field. For the crystal with size L y bigger 
than the radiation wavelength, k^L y 3> 1, the result is 
quite different: 



P 



uL y Ll\E z {Q,Q){ 
16tt 



[1 - Jo {k u L x )\ . 



(B6) 



Now the waves coming from opposite sides of the crystal 
along the y axis do not interfere with each other and radi- 
ation power becomes proportional to L y . For k u L x <C 1, 
we obtain the power proportional to fc 2 L 2 due to destruc- 
tive interference of the waves coming from opposite sides 
of the crystal along the x axis. If we put highly conduc- 
tive metallic screens separating the spaces x > L x and 
x < so that the edge x = radiates only into x < half- 
space, while that at x = L x radiates only into x > L x 
half-space (see Fig. 2), the interference will be eliminated. 



Such screens also double the radiation coming from one 
side. This can be demonstrated in the simplest way us- 
ing image technique^* radiation from the real electric 
currents induced at the screens is equivalent to radia- 
tion from the image magnetic current placed next to the 
original magnetic current. This leads to doubling of the 
effective magnetic current, M s — > 2M S , and quadruples 
the radiated power density. As the radiation now is lim- 
ited only by half-space, the total radiated power dou- 
bles. Therefore, in the presence of screens, the factor 
[1 — Jo(kujL x )] in Eq. (B6l has to be replaced by the fac- 



tor 2. This means that the screens strongly enhance the 
radiation induced by the homogeneous mode in the case 
k u L x <C 1. Such a design with screens for a crystal thin 
along the x axis was proposed in Ref. [23] This design 
gives the possibility of frequency tuning. In addition, 
heating is reduced due to small L x . However, the crystal 
should have a large number of layers to synchronize os- 
cillations in all junctions and work in the super-radiation 
regime. 

Next, we consider the fundamental cavity mode 
(m,n) = (1,0), more relevant for the subject of this pa- 
per. In this case, we obtain 



E _ L z exp^fc^r) 
£ 2 (1,0) ~ 4tit 



k kc k \ 

Px (n/L x ) 2 - fc| +Pjy fc^ J ^ ~ ex PH fc « L *)] i 1 + ex P i~ ik v L v) 



The components of the faraway electric field are given by the expressions 
Eg L z exp(ik^r) sin 2 9 — (tt/ci x ) 2 cos</ 



E<f, _ L z exp(zfc^r) (tt/Qx) 2 



47rr sin 2 #cos 2 <j> ~ (n/a x ) 2 sin sin 

cos 9 



1 + exp (— ia x sin 9 cos (f>)] [1 — exp {—ia y sin 9 sin <j))\ 



E Z (1,0) 47rr sin 2 (9 cos 2 - (ir/a x ) 2 sin 9 

with a x — k u L x and a y = k^Ly. The radiated power (B3l can be represented as 



- — — [1 + exp (—ia x sin 9 cos <p)\ [1 — exp (—ia y sin 9 sin <j))\ 



with 



cL 2 z \E z (l,0)\\ T , T , 
P = Ii,o{k^L x , k^Ly) 



7r ' 2 l" 7 ^ 2 [1 + cos (a x sin 9 cos <$)] [1 — cos (a y sin 9 sin (j>)] 
' d9 ^ 



(sin 2 #- (tt / a x ) 2 ) 2 cos 2 cf) + (tt / a x ) 4 cos 2 9 sin 2 <j) 
(sin 2 9 cos 2 4> — (Tr/a x ) 2 ) 2 



(B7) 



(B8) 



In the regime k^Ly 3> 1 this gives the following result Now we have a constructive interference of waves coming 



P : 



uL y Ll\E z {lM' 
16tt 



[1 + MkM] . (B9) 
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from opposite sides of the crystal along the x axis because 
electric field on these sides has opposite signs generating 
the same-sign magnetic fields. For such mode, screens do 
not influence much the radiation in the limit k u L x <C 1. 
The reduced parameter of radiation damping v r intro- 
duced in Eqs. (20 1 and (211 is related to the radiation 
power as 



P = u r L x L y L z £ -^\E z (l, 0)| 



16tt 



(BIO) 



The above results can also be straightforwardly gener- 
alized to the case when a stack is bounded by large-size 
ground plate at the bottom, 2 = 0. This is the case, for 
example, for the mesa fabricated on the top of bulk crys- 
tal. If we treat the ground plate as an ideal conductor, its 
influence can again be taken into account by the image 
technique^. This just leads to the doubling of the effec- 
tive magnetic current, M s — > 2M S , and to the doubling 
of the total radiated power P — > 2P. 
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